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ABSTRACT 

We have investigated the effects of uniform rotation and a specific model 
for differential rotation on the pulsation frequencies of 10 Mq stellar models. 
Uniform rotation decreases the frequencies for all modes. Differential rotation 
does not appear to have a significant effect on the frequencies, except for the 
most extreme differentially rotating models. In all cases, the large and small 
separations show the effects of rotation at lower velocities than do the individual 
frequencies. Unfortunately, to a certain extent, differential rotation mimics the 
effects of more rapid rotation, and only the presence of some specific observed 
frequencies with well identified modes will be able to uniquely constrain the 
internal rotation of pulsating stars. 

Subject headings: stars: oscillations, stars: rotation 

1. Introduction 

Observationally detected stellar pulsation frequencies can be used to place constraints 
on stellar models, giving us an improved understanding of the interior structure and evo- 
lution of stars. The most successful application has been the Sun, where the large num- 
ber of observed modes have placed strict constraints on parameters such as the helium 
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(IChristensen-Dalsgaard al. 



198S 



abundance fYl (jBasu &: Antial 12004 : lAntia fc Basul l2006l ) . the depth of the convection zone 



1991} 



an d the interior angular momentum distribu- 
tion (lEff-Darwich et al\ l2002l : ICouvidat et al\ l2003l ). Observations of pulsation frequen- 



cies of other stars continue to impro ve, particularly thr ough dedi cated satellites su ch as 
WIRE jHacking et a/.lll999h. M OST JWalker et a/.lboosh. CoRoT JSadin et a/.ll200lh and 
Kepler (jBasri. Borucki fc Koch 2005) as well as gro und-based networks such as STEPHI 
( Belmonte et al. 19931 ) and WET ( Nather et al. 1990 ). These improved observations, giving 
us long term coverage and improved accuracy, are the first steps in enabling other stars to 
be constrained in a similar manner to the Sun. Asteroseismology then, has the potential 
to answer a number of questions about the interior structure of stars throughout the HR 
diagram. 

One aspect of stellar structure which could be explored using asteroseismology is the 
internal rotation rate. It is theoretically possible for stars to rotate with angular velocity 
increasing or decreasing with distance from the rotation axis, and there is so me evidence tha t 
the latter may be true in massive main sequence stars, at least at the surface (jStoeckleylll968l ) . 
A third possibility is uniform rotation. It has been argued th at uniform ro tation is unrealistic 
based on observations of the Praesepe and Hyades clusters (ISmithlll97ll ). Of course, other, 
less well structured rotation laws are possible. However, there is little evidence in support 
of a specific rotation law, and the large uncertainties prevent any of the possibilities from 
being ruled out. We note in pas sing that the sola r rotation rate in the c onvection zone is 
primarily dependent on latitude (ISchou et a/.lll998uThompson et a/.ll2003l ) and thus cannot 
be described by a conservative rotation law. 



Recently, interferometric observations of Achernar (IDomiciano de Souza et al\ |2003| ) 
found that this star is far more oblate than is possible for a uniformly rotating star. This 
is true because uniformly rotating stars reach critical rotation before they have su fficient 
angular momentum to produce such an oblate object. However, iJackson et all (120051 ) noted 
that models in which the rotation rate increases inward from the surface can produce the 
oblateness observed for Achernar and still match the observed vsini. W hile further study ha s 
proposed that the oblateness may be due to a circumstellar envelope (ICarcofi et a/.ll2008l ). 
the original conclusion does raise the interesting question as to whether stars with rotation 
laws required to produce such an oblate shape exist, and if so, whether these rotation laws 
could be identified by possible pulsation modes. We investigate this possibility in this paper. 

Differential rotation with the rotation rate increasing inwards, as is considered in this 
paper, will have an impact on the deep interior structure of the star, provided the differential 
rotation is large enough. Rapid rotation in the outer layers of a star has little to no effect 
on the gravitational potential and core structure, as the envelope contains a relatively small 
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fraction of the stellar mass. In fact, many early attempts to model rotating stars assumed 
that the mass in the envelope was negligible and that the gravitational po tential in this 
region could be modelled using a Roche potential (ISackmann fc AnandlllQTOl ). However, it 
was recognized early on that this assumption was not always valid. Efforts to model a wider 
variety of rotati ng stars were made thr ough the implementation of the self-consistent field 
(SCF) method (lOstriker fc Marklll968r). or thr ough direct, 2D finite difference solutions to 



Poisson's equation (IClementI Il974j . 1 19781 . Il979l ). These methods allowed stars to be mod- 
elled with differential rotation, at least under certain circumstances. Concentrating angular 
momentum in the center, unlike uniform rotation, can produce enough distortion to affect 
the core, and consequently the evolution of the star. Only in this case can we produce 
a model with interior prop erties significantly different from the uniformly rotating model 
( ISackmann fc Anandl 119701 ). Even restricting ourselves to this type of differential rotation 
do es not narro w the possibilities significantly. The rotation could be shellular, as proposed 
by lZahnI (119921 ). or cylindrical (conservative rotation laws). In this paper, we have focused on 
conservative rotation laws, either with uniform rotation or with the rotation rate increasing 
towards the center of the star. Further discussion of our models can be found in fj2l In [|3]we 
consider the eigenfrequencies of rotating models as well as the large and small separations 
in ^and J|5] respectively. Our conclusions are summarized in ^ 



2. Numerical Models 



The st ellar models are computed using the 2D stellar structure code RO TORC (IDeupree 



1990l . Il995l). The code uses the OPAL opacities (llglesias fc Rogers! Il996l ) and equation of 
state (IRogers. Swenson fc Iglesiaslll996l ). Here we consider only 10 Mq ZAMS models with 
X=0.7, Z=0.02. These models solve the conservation equations of mass, momentum, energy, 
and hydrogen abundance along with Poissons equation for the gravitational potential on a 
two dimensional finite difference grid with the fractional surface equatorial radius and the 
colatitude as the independent variables. The surface equatorial radius is determined by 
requiring that the integral of the density over the volume of the model equals the stellar 
mass. The ZAMS models are taken to be time independent and static, except for the 
imposed rotation law, so that the mass, azimuthal momentum, and hydrogen composition 
conservation equations drop out. 

The only change required in the stellar evolution code for nonuniform cylindrical rotation 
laws is the addition of an extra term in the total potential (e.g., Tassoul 2000): 



(1) 
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= $ h / w^n{w) — i—r^dw' 

2 Jq dw' 

where VL is the rotation velocity (in radians per second) and vj is the distance from the 
rotation axis {xsinO, where x is the fractional surface equatorial radius and 6 is the colati- 
tude). The extra term is the last term on the right hand side of the equation. Performing 
the integral is straightforward for an analytically imposed rotation rate distribution in the 
ZAMS models. The total potential is used only to determine the surface location at each 
latitude by taking the surface to be an equipotential. Although defining a total potential 
requires a conservative rotation law, this is the only way in which a conservative rotation 
law is utilized in the stellar structure code. 

We have constructed uniformly rotating ZAMS models with rotation velocities between 
and 360 km s~^, with an approximate spacing of 30 km s~^. We have also computed 
a number of differentially rotating models at two values of the surface equa torial rotation 



velocity, 120 and 240 km s ^. The differential rotation law is as given by iJackson et al 



(120051): 



where /? is a parameter ranging from (uniform rotation) to 2, the maximum allowed for 
stability. The parameters a and Qo are used to impose the desired surface equatorial velocity 
and shape of the rotation law at small distances from the rotation axis. We have arbitrarily 
chosen a = 2. Figure [T] shows the rotation rate as a function of distance perpendicular to the 
rotation axis for a surface equatorial rotation velocity of 120 km s~^ and a surface equatorial 
radius for a uniformly rotating model at that speed. Increasing (3 increases the rotation rate 
close to the the rotation axis, including in the core of the star. Increasing angular momentum 
increases structural changes, and thus the structural changes increase with increasing /3. It 
is expected that increasing the rotation rate through increasing jS may in some ways mimic 
more rapid uniform rotation. 

One major result produced by significant rotation is an appreciable distortion of the 
surface of the model. We present the surface shape for a set of uniformly rotating models 
with surface equatorial velocities ranging from to 360 km s~^ in Figure [2l For each model 
the equatorial radius is taken to be unity. The ratio between the polar and equatorial radius 
decreases with increasing rotation, as the polar radius decreases slightly while the equatorial 
radius increases considerably. Differential rotation in which the rotation rate increases with 
decreasing distance from the rotation axis amplifies this effect. We present the surface shape 
for the differentially rotating models in Figure [31 The solid curves are for a surface equatorial 
velocity of 120 km s~^, while the dashed curves denote a surface equatorial velocity of 240 
km s~^. As the parameter j3 in Equation [2] increases, the fractional polar radius decreases. 
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The change in fractional radius with 13 is greatest at the pole and decreases towards the 
equator. Note that the fractional polar radius for a model rotating with a surface equatorial 
velocity of 120 km s~^ and a value of (3 of 1.8 has nearly the same fractional polar radius as 
a model uniformly rotating at 240 km s~^. 

We have increased the radial resolution of the static models by more than a factor of 
two over that used by Lovekin and Deupree (2008). The intent is to reduce the scatter 
and uncertainty in the pulsation mode calculations, particularly for the large and small 
separations. By and large this has been successful. 

The determination of the pulsation al proper t ies of these models is made using the linear 



adiabatic pulsation code developed by IClementI (Il998l ). We restrict our attention to input 



models with conservative rotation laws so that we can write the effective gravity [g) as the 
derivative of the total potential. The input models are axisymmetric spheroids, which allows 
us to assume a e*[^*+"^'^l time and azimuthal dependence. The equations to be linearized 
are the three components of the momentum conservation equation, the mass conservation 
equation, the adiabatic relation between the density (p) and the pressure (P), and Poisson's 
equation. The dependent variables are the three components of the linearized displacements 
(^), the linearized Eulerian displacements of the density pressure (5P), and the gravi- 
tational potential (50). We start with the perturbed momentum equation: 

where v is the velocity, which includes both the pulsational perturbations and the velocity 
distribution of the unperturbed model Here we assume this later velocity results exclusively 
from rotation, and is given by 

V = v(j) = Qrsin{6)(j). (4) 

Here Q is the rotation rate, and our assumption of a conservative rotation law requires it 
to be a function only of the distance from the rotation axis, w = rsin{6). We can write an 
expression for the Eulerian perturbation of the velocity in terms of the displacements because 
the Lagrangian time derivative of the displacement is the Lagrangian velocity perturbation. 
Thus, 

Sv = ^ + iv-y)^-i^-y)v. (5) 

Substituting this expressio n into the above perturbed momentum equation and proceed- 



ing along the same lines as IClementI (119981 ) leads to the following expression for the three 



individual components of the momentum equation: 



dr r do 



d6p 
dr 



-Agr{V-0 (6) 
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cr^^e + 2iail^ipCOs9 — wcosO 



a 



^ dr ' r do 



1 dSp 
r do 



(7) 



im 



■5p = 0, 



rsinO 

where a = uj+mQ, 6p = 6P/ p—S(j), A = c^dlnp/ — 1 and = TiP/ p. The quantity c is the 
adiabatic sound speed. The other equations required are as hsted in IClementI (1l998l ): mass 
conservation, the adiabatic relationship between the density and pressure perturbations, and 
Poisson's equation for the perturbed gravitational potential: 



5p = -V • pi 
5P = -ViPv -i-i-vp 

7^50 = -ATiG5p. 



(9) 
(10) 
(11) 



Differential rotation has introduced two new terms into each of the radial and latitu- 
dinal momentum equations. The other change is that cr, the eigenfrequency in the local 
rotating frame, is no longer constant except for axisymmetric modes. This does not alter 
the solution algorithm, however, because Clement's (1998) approach was to select a (which 
now becomes selecting u and calculating a locally as needed), solve all the equations and 
see if a discriminant was satisfied. Neither the added terms nor the nature of the eigenvalue 
requires any further mod ification t o the approach, and the derivation of the final equations 
proceeds as described by IClementI (119981 ). 



Using these equations we can calculate the pulsation properties of the stellar models 
on a 2D finite difference grid. This is done through a change of variables, factoring out the 
behavior of ^e, and 50 near the boundaries to eliminate singularities. The coefficients of 
these equations can be put in a band-diagonal matrix and solved. NRO can include up to nine 
angular zones in the eigenfunction solution. This gives the solution at N angles, where N is 
the number of angular zones, which can subsequently be decomposed into the contributions 
of individual spherical harmonics through the use of Fouri er transforms. Thro u ghout this 
paper, we have used N = 6. Based on the calculations of iLovekin fc Deupred (120081 ) . six 
spherical harmonics is sufficient to accurately calculate the eigenfrequencies for the most 
rapidly rotating models discussed here. Indeed, we have performed a few test calculations 
with N = 8 and have found that the effect on the frequencies is small, typically a few 
hundredths of a percent. 



As discussed in iLovekin fc Deupred (120081 ). NRO, combined with stellar structure mod- 
els from ROTORC, allows us to calculate the pulsation frequencies for rotating stars without 
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making any a priori assumptions about the structure, except that the rotation law is con- 
servative for NR O. For further discussion of t he method of solution used by NRO, refer to 
Clement Jl998h 



or 



Lovekin fc Deupred ((20081). 



With spherical stellar models, the radial and angular components of the perturbations 
separate, and the angular part can be expressed as a spherical harmonic with specific values 
of the quantum numbers, / and m. 

For rotating stars, the eigenfunction solution is not a single spherical harmonic, and / 
is not a valid quantum number. Indeed, NRO uses I only to specify the parity of the mode 
being calculated, and includes the first k even or odd spherical harmonics, where k is the 
number of angles included. We identify modes using a quantum number lo, which is the / of 
the mode in the non-rotating model to which a given mode can be traced. For spheroids, m 
remains a valid quantum number. As in iLovekin fc Deupred (|2008[ ) , we restrict ourselves to 
axisymmetric modes (m = 0) and modes with small radial quantum number (n). 



3. Relative frequencies 

In this paper we consider low order axisymmetric modes for = 0, 1, 2 and 3. These 
modes are expected to have the highest amplitudes and the smallest cancellation effects 
across the visible surface of the star, and are hence expected to be the most easily visible. 
Our structural models cover velocities from to 360 km s~^ and for two velocities, 120 and 
240 km s~^, we have calculated differentially rotating models with /3 varying from to 2.0. 
Tracing the individual modes becomes very difficult above rotation velocities of 360 km s~^ 
and for some higher values of f3, and this represents a practical limit to our study. Although 
the frequencies can be calculated at these velocities, the resulting eigenf unctions are a mix of 
six spherical harmonics, and no single harmonic dominates. As it is very difficult to reliably 
assign a value of lo to these modes, we exclude them from our analysis. It is probably feasible 
to trace the modes accurately, but this could require an extremely fine rotational velocity 
grid (1-5 km s~^). We decided not to pursue this for this exploratory work. For differentially 
rotating models, the limits are (3 = 1.8 for the 120 km s~^ model and (3 = 1.0 for the 240 km 
s~^ model. Based on the curves shown in Figure [1], it appears that the limit is related to the 
angular velocity near the rotation axis. The curve representing /3 = 1.0 has approximately 
half the value at the center of the f3 = 2.0 curve. Therefore, if we double the velocities 
everywhere, the limiting f3 should move from 1.8 at 120 km s~^ to 1.0 at 240 km s~^, which 
corresponds to approximately the same angular velocity near the rotation axis. 
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3.1. Uniform Rotation 



The trends produced by tracing a given mode tlirough increases in rotation velocity are 
illustrated in Figure H] for the = 2 mode, which shows the eigenfrequencies normahzed by 
the non-rotating frequency for each mode. Overall, the trends we find for frequency agree 
with those calculated by previous work (iLignieres et al\ l2006l ) . These authors find that 
the frequencies decrease as one increases the rotat ion rate, with highe r frequ ency modes 
decreasing more than lower modes. As discussed in iLovekin fc Deupred (120081 ). our results 
at low to moderate rotation rates are also consistent with the frequency trends predicted by 
second order perturbation theory (see for example, Saio 1981). 

We have increased the radial resolution of the outer 30% of the radius of the static 



models by more than a factor of two over that used by iLovekin &: Deupred ( 120081 ). This 
produces a radial zoning finer than that currently allowed by the pulsation code, so further 
increases in radial resolution in the 2D structure models will only be effective if the pulsation 
code is modified to allow more radial zones. The intent of the r nodified zoning is to reduc e 
the scatter and uncertainty in the mode calculations evident in ILovekin fc Deupred (120081 ). 
Figure H] shows that a reasonable estimate of our accuracy for the frequencies is a very 
few tenths of a percent, although there are still a few frequencies, most commonly but not 
exclusively for the higher radial orders and higher rotation rates, which do not fit within this 
limit. One might expect the radial resolution near the rotation axis and at mid latitudes to be 
less than for lower rotation rates because the fractional radius at these latitudes compared to 
the equator is lower. The accuracy of the small separation appears well within a /xHz, while 
the large separation does show variations on the order of one /iHz, particularly at higher 
rotation velocities and for higher radial order modes. This is compatible with the notion that 
the radial resolution near the surface could continue to benefit from refinement. However, 
these uncertainties do not disguise trends in the results, even in the large separation, with 
respect to rotation rate or the rotation law, and we consider these trends significant. 

One interesting line of inquiry is whether there is some analog to the period - mean 
density relation which allows interpolation of eigenfrequencies as functions of models and 
rotation rates. Specifically, we have examined if there is a physically meaningful radius 
which can be used in the period - mean density relation 



Q 




(12) 



(where M and R are in solar units, and P is the period) that would allow Q to be approx- 
imately constant as a function of the rotation rate. The comparatively small changes in 
the eigenfrequencies shown in Figure H] suggest that the surface equatorial radius, with its 
fairly rapid increase as a function of rotation, will not keep Q approximately constant, and it 
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does not. The same is true for an average radius, defined as eitlier a straight average or the 
effective radius required to contain the total volume of the model. The polar radius would 
be more promising because it only slowly varies with the rotation rate, but it actually de- 
creases slightly as the rotation rate increases. This is the wrong direction to keep Q constant 
because the frequency decreases as well. Because the polar radius decreases slightly and 
the equatorial radius increases appreciably with increasing rotation, one might guess there 
would be some latitude at which the radius increases at a rate that nearly offsets the rate 
of period increase. This is true, and occurs at a colatitude of 40 degrees. It is not obvious 
that this has any physical significance because it is difficult to associate any specific meaning 
to the radius at this latitude. We present the pulsation constant for two definitions of an 
effective radius in Figure [51 One way is to use the radius of a sphere with the same volume 
as the model. The other uses the radius at a colatitude of 40°. For comparison we show a 
"pulsation constant" that would exist if we used the mean density of the nonrotating model. 
We include this to give an idea of the size of the effect. Interestingly, the variation in the 
pulsation constant is significantly larger allowing the mean density to be determined by the 
total volume of the rotating model (the mass is the same for all models) than it is when the 
mean density is assumed to be that of the nonrotating model. 

The frequencies can also be changed by the mass or evolutionary state of the star, 
producing trends that could potentially be confused with rotational effects. We wish to 
determine how closely the frequencies of a rotating model can be mimicked by a non-rotating 
model. First we calculated the Q values for each model in the 10 and 12 Mq non-rotating 
models. For each mode, we then took the mean of the Q of the two models. We used this 
average Q for the radial fundamental mode and the frequency of the radial fundamental 
mode for the model rotating at 150 km s~^ to calculate a mean density. This corresponds to 
the mean density of a non-rotating model of unknown mass and radius pulsating in the radial 
fundamental mode with the same frequency as the 10 Mq model rotating at 150 km s~^. The 
mean density found this way and the average Q's for the other Zq's can be used to predict 
the other frequencies of this presumed non-rotating model. When these frequencies are 
compared to the calculated frequencies for the rotating model, the differences are significant. 
Using Q to calculate the frequencies in this way forces the radial fundamental mode to 
have the same frequency, so the differences between frequencies should be solely a result of 
rotation. The frequencies predicted for the lo = and 2 modes are larger by 1-5%, with the 
differences increasing for higher order modes. At the same time, the frequencies predicted 
for the lo = 1 and 3 modes are smaller, by as much as 15 % for the lo =1 pi mode. As 
the radial order increases, the differences between the rotating model and the non-rotating 
calculation decreases for the /o = 1 modes, but increases for the /o = 3 modes. The size and 
direction of these trends implies that the pulsation spectrum of a rotating model is unlikely 
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to be confused with the pulsation spectrum of a more massive non-rotating model. It also 
suggests that rotation must be included in the calculations if observations indicate it might 
be present even at this moderate amount. 

We have also evolved a single non-rotating model, and compared the ZAMS model with 
one part way through the main sequence evolution (Xc = 0.47). In this case, the frequencies 
decreased sufficiently, even for a model with a large remaining core hydrogen fraction, that 
confusion seems unlikely. 

3.2. Differential Rotation 

We have studied the change in the frequencies of IOMq ZAMS models differentially 
rotating at 120 km s~^ and 240 km s~^. Overall, the frequencies increase for lo = and 1, 
and decrease for lo = 2 and 3, a trend seen at both 120 km s~^ and 240 km s~^. Our results 
for lo =0 are shown in Figure M for the fundamental, IH, 2H and 3H modes for a model 
rotating at 120 km s~^. In this case, the frequency changes are largest for the 3H modes, 
but are noticeable for all modes by /3 ~ 1. Similar trends are found for the other values of 
lo considered here. Still, the differences remain relatively small, and it seems unlikely that 
even extreme differential rotation with this surface rotation velocity will be detectable using 
the values of the eigenfrequencies alone. 

The frequency results for the 240 km s~^model, shown in Figure [7] for lo = 0, are slightly 
more promising. Although we were unable to reliably identify modes above (3 = 1.0, the 
frequencies already differ by more than 1% hj f3 = 1.0 for the IH mode, and it seems the 
differences would be noticeable hj (3 = 0.4. If this trend continues, as seems likely at least 
for the F and IH modes, the frequency differences should be large enough to be detectable 
in these more rapidly rotating stars. As noted above, for some modes differential rotation 
causes the frequencies to increase as in Figures [7] and [H while for others the frequencies 
decrease, as in Figure M These plots do not include the 3H mode, as we found that the 
scatter in this mode remained a significant fraction of the variation, despite the improved 
radial zoning, and so have chosen not to include it in our discussion. 

The effects of differential rotation compared with uniform rotation are shown in Figure 
in] for the lo = 2 p2 mode. Differential rotation can change the frequencies by about a percent 
above and beyond the difference predicted based on surface equatorial velocity alone. The 
differences are small; about 1% for the most extreme differentially rotating model at 120 
km s^^. Based on the frequencies we have calculated, it may be possible to discriminate 
between uniform and this type of differential rotation, given the right combination of properly 
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identified frequencies. Since frequencies increase relative to tlie uniformly rotating case for 
some lo, and decrease for others, these differences could be used to constrain the rotation. 
This would require a star with a few positively identified modes, some of which were either 
= or 1, and some of which were either = 2 or 3. The number of modes required and 
the challenges presented by accurate mode identification in massive main sequence stars may 
make this extremely difficult in practice. 



4. Large Separations 

If one considers two stellar models that are in the same evolutionary phase, and appear 
reasonably close to each other in the HR diagram, the frequencies can be approximately 
determined from the relevant pulsation constant, Q. As the mass and radius change with 
position in the HR diagram, so will the frequencies. It is expected that for small changes 
in mass and radius the frequency differences (either frequency separations or ratios) will 
change, like Q, much more slowly than the individual frequencies. As a result, the large 
and small separations are probably more useful than individual frequencies as they are less 
sensitive to small changes in the models. 



4.1. Uniform Rotation 

The large separation, defined as 

Az/, = Ui^n+l - ^l,n (13) 

can provide information about the outer layers of the stellar envelope. The large separation 
for the lo = mode is shown as a function of the rotation rate in Figure [TOl We note 
that the overall trend of the large separation for this mode is to decrease as the rotation 
rate increases. At other lo the trend is the same and the large separation decreases for 
every pair of modes considered. The magnitude of the decrease in large separation does 
increase slightly with increasing lo, as can be seen by comparing Figures [10] and [TTl For 
nonrotating ZAMS models, both the frequencies and the large separations decrease as the 
mass (hence the radius, luminosity and effective temperature) increases. The decrease in 
the large separation occurs not only because the frequencies decrease, but also because the 
period ratios increase for increasing ZAMS mass. However, for stars observed approximately 
equator on, rotation decreases the perceived luminosity and effective temperature. This 
offset between the perceived luminosity and temperature and the large separation may be 
useful as a rotation discriminant. Of course, stars observed nearly pole on show an increase 
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in perceived luminosity and effective temperature, which is in hne with the decreasing large 
separation as the rotation rate increases. Any discriminant of rotation may only be a matter 
of degree for low inclination objects. 

We noted in section 13.11 that the decrease in the frequencies with increasing rotation 
cannot be explained purely by the decrease in the mean density (i.e., a constant Q). The mean 
density decreases faster than the pulsation periods increase as the rotation rate increases. 
Interestingly, the decrease in the large separation for /q = in Figure [TD] is almost e ntirely 
offset by t he mean deri s ity so that Au{pq/ pY^"^ is nearly constant, as discussed by lUlrich 
( 119861 ) and lReese et al\ (120081 ). The mean density does not offset the steeper decline in the 
large separation for the /o = 2 modes shown in Figure [TTJ 



4.2. Differential Rotation 

The large separations provide information about the region near the surface of the star, 
while the frequencies provide information on more global properties. As discussed above, the 
separations are less sensitive to small changes in the mass or radius of the star, but since they 
probe the surface region, may provide information about changes in this region as a result 
of rotation. A comparison of Figures [2] and [3] shows that the polar radius is significantly 
more affected by differential rotation than the radius at lower latitudes. This kind of effect 
may be detectable using the large separations. Indeed, based on the results in the previous 
section, we expect there to be significant changes in the period differences, as we have found 
that differential rotation can sometimes introduce a significant shift in only one or two of 
the harmonics. 

Although the individual frequencies can show significant relative differences, this does 
not carry over to the large separations. The large separation for the lo = modes of a dif- 
ferentially rotating model with surface equatorial velocity of 120 km s~^are shown in Figure 
[121 The large separations shown in this plot show very little change with increasing (3, much 
less than the differences shown in Figure [TUl The same lack of variation is seen for all modes. 
Given that the large separation probes the surface regions, this might be regarded as some- 
what surprising because changing the (3 does change the surface configuration, particularly 
near the rotation axis. 

At 240 km s~^, the large separations, shown in Figure [121 are again quite constant over 
the region shown. Most separations either remain constant or show a slight increase, at least 
to (3 = 0.6, at which point some of the higher order separations decrease slightly. Again, this 
is different from the trend seen in the uniformly rotating models. Particularly at high f3, the 



- 13 - 



separations moving in different directions may allow constraints to be placed on observed 
stars. The differences begin to become noticeable at (3 ~ 0.6 for most of the high order 
modes considered. However, for most modes the large separation never differs by more than 
a few //Hz. As for the frequencies, it seems that the large separations are unlikely to produce 
any very refined constraints on the internal rotation rate, at least for this particular rotation 
law. 



5. Small Separation 



Asymptotic theory (lTassoullll980l ). which predicts that the large separation should be 
approximately constant as n gets large, also predicts near degeneracy between modes with 
the same value of n + 1/2: 

T^l,n - l^l+2,n-l (14) 

The deviations from this degeneracy are defined as the small separation: 



di 



-(4/ + 6) 



dc dr 



4vr2zy, „ Jq dr r 



(15) 



where c is the sound speed. The sound speed changes in the core are sufficiently large that 
the (1/r) variation dominates the integral on the right hand side of Equation [151 s-iid hence 
the small separation is dominated by the structure in the core. 



5.1. Uniform Rotation 

At slow uniform rotation, the size and shape of the convective core is nearly unaltered by 
the rotation, and one would expect the effects on the small separation to be minimal. Figure 
m shows this to be true, but also shows that the small separation for higher n increases 
markedly with the rotation rate once the rotation exceeds approximately 150 km s^^. There 
are slight changes to both the shape and relative size of the convective core with rotation, 
although the absolute mass and radius of the core change only slightly. It is not obvious why 



the sm all separation increases so markedly. This effect has also been noted by lLignieres et al. 



mm 



Small separations are frequently used as probes of th e core structure of stars, and can 



be used to constrain overshooting and core composition (ISoriano fc Vauclairl l2008l ) . Their 
results indicate that convective core overshooting causes a slight decrease in the small sep- 
arations, with the effect becoming more pronounced as the star evolves. This slight trend 
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is opposite to that produced by at least moderate rotation, which appreciably increases the 
small separation. Clearly, this is a situation in which caution must be exercised when using 
observed modes to constrain conditions deep in the stellar interior. 

5.2. Differential Rotation 

For differentially rotating models, the overall trend is the same, with small separations 
increasing with increasing differential rotation. However, Figures [15] and [16] show that the 
variation in small separation is much less than for uniformly rotating models. The trends 
are consistent with the relationship between the effects of (3 and those of increasing the 
uniform rotation rate as shown in Figure [9] The effects of increasing (3 on the convective 
core mimic to some extent those of increasing the uniform rotation rate, although high values 
of j3 do make the convective core more oblate. The different effects of j3 on the large and 
small separations is understandable in that increasing (3 increases the rotation and its effects 
near the rotation axis, and this certainly includes the convective core. However, this trend 
of increasing (3 producing similar trends to increasing uniform rotation rate does not give 
us confidence that we have a useful tool for diagnosing a rotation law of the kind we have 
considered through the small separation. 

6. Conclusion 

We have investigated the effects of uniform and differential rotation on pulsational eigen- 
frequencies. For uniformly rotating models, we have found that the frequencies decrease as 
the rotation rate increases for all values of lo and n considered here, although the rate of de- 
crease varies with the mode in question. While this frequency behavior is expected assuming 
the period-mean density relation applies, the frequency changes are much smaller than the 
period- mean density would suggest. We do find a pulsation constant being approximately 
constant if we use the surface radius at a colatitude of about 40° in the period-mean density 
relation, although this radius does not represent the mean density. 

For the differential rotation law considered here, we find the frequencies at a given ve- 
locity may either increase or decrease, depending on Zq, with increasing differential rotation, 
relative to the uniformly rotating model. However, the overall effects in all cases are com- 
paratively small, with maximum differences typically on the order of 1% when compared to 
the uniformly rotating case. 

Uniform rotation decreases the large separation by several /iHz (< 10) over the entire 
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range of rotation (0 - 360 km s~^) considered here. The large separation was virtually 
unchanged (< 1 /xHz) from that of uniform rotation for the range of differential rotation 
parameters considered here, despite the noticeable change in the surface shape. Although 
this change in shape is noticeable, it is still considerably smaller than the change produced 
by uniform rotation. Both uniform and differential rotation increase the small separation. 
The small separation can change markedly over the range of uniform rotation considered, 
while the dependence of the small separation on the rotation profile is more modest but 
not inconsistent with the other effects produced when comparing uniform and differential 
rotation. The effects of rotation on the frequencies and separations are generally large enough 
that rotation must be considered in the asteroseismology of these upper main sequence stars. 
While the precise rotation rate at which one must be concerned with rotation depends on 
the level of accuracy achievable, it is certainly no larger than 100 km s~^ for our 10 Mq 
ZAMS models. 

Although we have shown there can be significant differences in the pulsation properties 
of rotating stars, it is not clear that these results can actually be used to constrain the 
interior rotation rate. Given the possible combinations of effects from the rotation rate and 
distribution, the mass of the star, convective overshoot, evolutionary stage, etc, it seems 
unlikely that pulsation properties will give a unique solution, particularly if the number of 
observed modes is modest or cannot be properly identified. However, we have found that 
some combinations of modes constrain some of these properties. 

This research was supported by the National Science and Engineering Research Council 
of Canada through a Discovery grant and a graduate fellowship. Computational facilities 
were provided with grants from the Canadian Foundation for Innovation and the Nova Scotia 
Innovation Research Trust. 
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Fig. 1. — Rotation law used in differentially rotating models (Eqn. [2]). Curves show from 
bottom to top the rotation law for /3 = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8 and 2.0 for a 
model with surface equatorial velocity of 120 km s~^. 
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Fig. 2. — Surface shape for uniformly rotating models. The polar radius decreases relative 
to the equatorial radius as rotation increases from km s~^ to 360 km s~^. 
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Fig. 3. — Surface shape for differentially rotating models at 120 km s~^ (solid) and 240 km 
(dashed). As the rotation rate close to the rotation axis increases (increasing /3), the 
polar radius decreases relative to the uniformly rotating case. 
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Fig. 4. — The first four harmonics of the lo — mode for a uniformly rotating 10 M0 
model as a function of the rotation rate. The four curves represent the frequencies for the / 
(diamond), pi (X), p2 (square) and p3 (triangle) modes. 
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Fig. 5. — Plot of the pulsation constant for the radius of a sphere with the same volume 
as the rotating model (bottom) and the radius at colatitude of 40° (middle) . Also shown is 
the pulsation constant assuming there is no change to the mean density as the rotation rate 
increases (dashed curve). 
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Fig. 6. — The fundamental and first three overtones of the lo = mode for a model rotating 
at 120 km s^^as a function of the differential rotation parameter (3 (see Equation [2]). The 
four curves represent the frequencies for the fundamental (diamond), IH (x), 2H (square) 
and 3H (triangle) modes. 
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Fig. 7. — Relative frequencies of the = modes versus differential rotation parameter 
13 for a model rotating at 240 km s^^. The curves show the relative frequencies for the 
fundamental (diamond), IH (x) and 2H (square) modes. 
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Fig. 8. — Relative frequencies of the lo — i modes versus differential rotation parameter for 
a model rotating at 240 km s~^. Shown are the pi (diamond), p2 (x) and (square) modes. 
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Fig. 9. — The relative effects of differential rotation for the Iq — 2 p2 mode. The frequencies 
for differentially rotating models as a function of (3 at 120 km s~^and 240 km s~^ (squares) 
are superimposed on the uniformly rotating frequencies (diamonds). 
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Fig. 10. — The large separation between the 3H and 2H (square), the 2H and IH (x), and 
IH and F modes (diamonds) for modes with = 0. 
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Fig. 11. — The large separation for the lo — "2 modes as a function of rotation velocity. 
Shown are the separations between the pi and / modes (diamonds), the p2 and pi modes 
(x) and the ps and p2 modes (squares). 
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Fig. 12. — The large separation for the lo = modes of a differentially rotating model 
with surface equatorial velocity of 120 km s~^, plotted as a function of differential rotation 
parameter /3. Symbols are the same as Figure fTOl 
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Fig. 13. — The large separation of the /o = modes for a differentially rotating model with 
surface equatorial velocity of 240 km s~^. Symbols are the same as Figure [TOl 
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Fig. 14. — Small separation for the lo — and 2 modes as a function of surface equatorial 
velocity. Shown are the separations between the lo — 0, 3H - lo — 2, p2 modes (squares), lo 
— 0, 2R - lo — 2, pi modes (x) and lo = 0, IH - = 2, / modes (diamonds). 
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Fig. 15. — Small separations for the even modes for a differentially rotating model with sur- 
face equatorial velocity of 120 km s~^, plotted as a function of differential rotation parameter 
p. Symbols are defined as in Figure [TH 
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